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Abstract 


The Downingia yina species complex (Campanvilaceae), centered in northern California and 
southern Oregon, currently contains three morphologically distinguished species: D. yina, D. elegans, 
and D. hacigalupii. This complex of species is notable for high levels of morphological and cytological 
variation, with chromosome counts of n = 6, 8, 10, and 12. Molecular evidence suggests three main 
clades within this complex, corresponding more with cytological variation than with morphological 
variation. Additionally, the molecular evidence suggests a phylogeographic pattern associated with the 
Cascade Ranges, where members of the clade characterized by chromosome counts of n = 6, 8, and 10 
are distributed primarily to the west of the Cascades while members of the clade characterized by 
chromosome counts of n = 12 are distributed primarily to the east. A third clade characterized hy n = 
10 is localized in the Lake of the Woods region of southern Oregon. Evidence from morphological, 
cytological, interfertility, and molecular data was used to re-examine the delimitation of species within 
this complex. Downingia elegans and D. bacigalnpii are maintained, while D. yina is split into three 
morphologically cryptic species (D. yina, D. willaniettensis, D. pulcherrinia) that do not form a clade. 
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The Downingia yina species complex is a 
monophyletic group (Schultheis 2001) comprising 
D. yina Applegate, D. hacigalupii Weiler, and D. 
elegans (Lindl.) Torr. The species complex 
represents a cytologically and morphologically 
variable group centered in northern California 
and southern Oregon. Chromosome numbers 
within the complex include n = 12 in Z). 

bacigalupii, n = 10 in D. elegans, and races of n 
= 6, 8, 10 and 12 in D. yina (Weiler 1962; Foster 
1972; Lammers 1993). Morphologically, both D. 
bacigalupii and D. elegans are distinguished from 
D. yina by an exserted staminal column with a 
sharp bend between the anthers and filament, and 
by the concave oval-shaped lower corolla lip with 
relatively parallel corolla lobes. Downingia baci¬ 
galupii can be distinguished from D. elegans by 
the corolla’s lighter shade of purple and by the 
yellow pigmentation in the corolla throat, a 
feature also found in D. yina. 

Morphological variation within D. yina has led 
some workers to recognize additional species or 
infraspecific taxa. Downingia yina was described 
by Applegate (1929) from a localized region of 
the southern Cascade Ranges in Klamath Co., 
Oregon. Shortly thereafter. Peck (1934, 1937) 
described two additional larger flowered species: 
D. willaniettensis Peck from the Willamette 
Valley of Oregon, and D. puleherrinia Peck from 
eastern Oregon. In the first monograph of the 
genus, McVaugh (1941) recognized D. yina and 
D. willaniettensis, including D. puleherrinia in the 
latter. McVaugh noted (1941), however, that D. 
yina and D. willaniettensis were not readily 
distinguishable, and ultimately treated them as 


varieties within D. yina, var. yina and var. major I 
McVaugh, respectively (McVaugh 1943). He j 
distinguished the two varieties based on fruit ' 
characteristics (fusiform with hyaline lines in var. 
yina; subulate without hyaline lines in var. major), 
plant stature (larger and more erect in var. 
major), and geographic location of the popula¬ 
tions. Weiler (1962) found that the differences 
described between D. yina and D. willaniettensis 
were not maintained under greenhouse condi¬ 
tions. He accordingly recognized only D. yina 
with no infraspecific taxa, although noting that ; 
fresh material of D. pulclierrima was not exam¬ 
ined. Weiler (1962) also noted that individuals of 
D. yina sensu lato tended to be decumbent to the 
west of the Cascade Ranges, and erect to the east. 
Foster (1972) was unable to find consistent 
morphological differences to correspond with 
cytological races in D. yina, but did note an 
ecological trend. She observed that D. yina 
chromosome races n = 6, 8, and 10 are found 
in habitats characterized by Kiichler (1964) as 
Oregon-oak woodland or cedar-hemlock-Doug- 
las fir mosaic while the D. yina chromosome race 
11 = 12 is found in California mixed evergreen i 
forest and juniper-steppe woodland, as charac¬ 
terized by Kiichler (1964). Both Foster (1972) 
and Ayers (1993) followed Weiler (1962) in 
recognizing only D. yina. 

The present study emerged largely from a 
systematic investigation of the genus Downingia, 
in which molecular data unexpectedly suggested 
the existence of morphologically cryptic lineages 
within D. yina (Schultheis 2001), corresponding 
in part to infraspecific taxa previously recog- ' 
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Fig. 1. Map of the northwestern USA showing localities of samples included in this study. Each symbol may 
represent one or multiple samples from the vicinity. The dashed line roughly corresponds to the geographic barrier 
created by the Cascade Range. Triangles = Donningia elegans. Squares = D. bacigalupii. Circles = samples 
previously included in D. yina. Filled circles = now assigned to D. wiUametteiisis. Open circles = now assigned to D. 
pulcherrima. Circles with line through center = now assigned to D. yina sensu strictii. Clade I, Clade II and Clade 
III refer to clades identified in phylogenetic analyses. 

nized. The situation was further eomplieated 
(Schultheis 2001) by the apparent para- or 
polyphyly of D. yina with respeet to D. elegans 
and D. bacigalupii. The D. yina speeies eomplex 
thus represents a mixture of morphologically 
cryptic and morphologically distinctive lineages 
that may not correspond to the species currently 
recognized (Ayers 1993). The aim of this study 
was to further investigate the relationships and 
circumscriptions of D. hacigalupii, D. elegans and 
D. yina using morphological data, additional 
nuclear and chloroplast molecular sequence data 
to supplement Schultheis (2001), and available 


cytological and interfertility data (Weiler 1962; 
Foster 1972). 

Methods 
Taxon Sampling 

Collections were made from throughout the 
range of D. elegans, D. bacigalupii, and D. yina 
(Appendix 1; Fig. 1). Herbarium collections 
provided important supplemental material. 
Downingia bicornuta A.Gray, D. concolor E. 
Greene, D. cuspidata (E. Greene) Rattan, D. 
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Table 1. Characters Used in Morphological Analyses of the Dowingia yin a Complex. Characters 1- i 
10 are quantitative and measured in millimeters, characters 11-14 are qualitative, and characters 15-19 are ratios. 

‘ Characters used in cladistic analyses, with character states noted in brackets. 


Character 

number 

Character 

Character definition and how assessed 

1. 

sepal 

dorsal sepal, length 

2. 

back slit 

corolla base to dorsal slit, length; equivalent to height of corolla tube along 
dorsal surface 

3. 

side slit 

corolla base to lateral slit, length; equivalent to height of corolla tube along 
lateral surface 

4. 

upper lobe 

upper corolla lobes, length 

5. 

lower lobe 

lower corolla lip, length 

6. 

filament' 

filament tube, length (<6 mm [0], >6 mm [1]) 

7. 

anther 

anther tube, length 

8. 

anther angle' 

angle between anther and filament tubes (<50 [0], >70 [1]) 

9. 

lower angle 

angle of divergence between lobes of the lower corolla lip 

10. 

horns' 

anther horns, length (<0.62 mm [0], >0.62 mm [1]); refers to triangular 
projections on each of the two smaller anthers 

11. 

anther back 

trichomes on anther dorsal surface: abundant (0), few (1), none (2) 

12. 

upper lobe orientation 

upper corolla lobes, orientation: parallel (0), intermediate (1), divergent (2) 

13. 

yellow' 

yellow on lower corolla lobe: present (0), absent (1) 

14. 

lower lobe shape 

lower corolla lobe, shape: acute (0), intermediate (1), mucronate (2) 

15. 

filament/anther 

filament length/Anther length 

16. 

sideslit/backslit' 

length of lateral slit/Length of dorsal slit (sO.6 [0], <0.6 [1]) 

17. 

filament/backslit' 

filament length/Length of dorsal slit (^1 [0], 1-2 [1], >2 [2]) 

18. 

upper/lower lobe 

upper lobes length/Lower lip length 

19. 

backslit/upper lobe 

length of dorsal slit/Upper lobes length 


luontana (E. Greene) Rattan, and D. ornatissima 
E. Greene were chosen as outgroup taxa based on 
previous phylogenetic analyses within the genus 
(Schultheis 2001). 

Molecular 

Generation of sequence data. Extraction of total 
DNA from 24 samples first reported in Schultheis 
(2001) and 9 new samples (Appendix 1) involved 
use of either the CTAB protocol of Doyle and 
Doyle (1987) or Hillis et al. (1996) with minor 
modifications (Schultheis 2001), or use of Qiagen 
DNeasy Plant mini kits following manufacturer’s 
instructions. Most plant tissue samples were 
stored in a cooler while in the field and transferred 
to a — 80C freezer within one week of collection. 
Voucher specimens were either the same plant 
from which tissue for DNA extraction was taken, 
or were other plants from the same site. 

Sequence data were generated from the nuclear 
18S—26S rDNA internal transcribed spacer (ITS) and 
the chloroplast 3'//7/K intron. Amplification and 
sequencing methods changed during the course of 
the project, as new techniques became available. 
Single-stranded DNAs of ITS 1 and ITS 2 were 
generated, purified, and manually sequenced 
following Baldwin (1992). Double-stranded 
DNAs of ITS 1, ITS 2 and the STrnK intron were 
generated, purified, and sequenced using automat¬ 
ed sequencing technology following Schultheis 
(2001). Sequences are deposited in Genbank. 


Sequence alignment. All alignments were visual. 
Sites coded with “?” or with an lUPAC-IUB 
ambiguity code represent basepairs where se¬ 
quence produced with neither primer produced a 
sufficiently strong or clear signal for confident 
basepair assignment. Indels, coded as were 
treated as missing data. Two regions were 
excluded from the ITS dataset due to ambiguous 
sequence alignment (positions 132-138 and 284- 
291 of the aligned ITS data set). 

Evaluation of sequence data. Separate and 
combined analyses using a parsimony criterion 
were conducted for ITS and y trn¥^ intron data. 
All analyses employed heuristic searches with 
10,000 replicates of random sequence addition 
and tree-bisection-reconnection (TBR) branch 
swapping. Conservative estimates of clade sup¬ 
port were assessed using 10,000 replicates of the 
“fast” bootstrap option in PAUP 4.0b5. Decay 
analyses (Donoghue et al. 1992; Bremer 1994) 
using Autodecay (Eriksson 1998) were conducted 
for the y trn¥^ intron and the combined molec¬ 
ular analyses. 

Morphology 

Fresh and/or herbarium material was exam¬ 
ined from 80 localities (Appendix 1, Fig. 1) and 
450 flowers. Nineteen characters were included 
for phenetic analyses, including 10 quantitative, 4 
qualitative, and 5 ratio characters (Table 1). All 
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characters are floral, because vegetative charac¬ 
ters are not generally useful for distinguishing 
species of Domiingia. Characters were observed 
or measured against a ruler under a dissecting 
scope, except for anther horn length which was 
measured with an ocular micrometer. 

Morphometric analyses. Analyses of variance 
were conducted to identify characters differing 
significantly among the three currently recognized 
species, and Tukey tests were used to identify 
which species differed. The same was done within 
D. yiiia for the three groups identified by 
molecular analyses (see results). For multivariate 
analyses, a data matrix was created containing the 
average value for each character from each 
collection locality. Multivariate analyses included 
cluster analysis using Euclidean distances and 
single linkage, discriminant function analysis, and 
Principal Components Analysis (PCA), the latter 
using a matrix standardized so that each character 
had a mean of zero and a standard deviation of 
one. All statistical analyses were performed with 
SYSTAT 5.2.1. 

Cladistic analyses. One qualitative and five 
quantitative characters (indicated in Table 1) were 
used in a cladistic analysis of the 26 populations 
for which molecular data were also available, plus 
one population per outgroup taxon. Phylogenetic 
analyses using a parsimony criterion were con¬ 
ducted with PA UP 3.1.1 (Swofford 1993) or 
PAUP *4.0b5. The analysis employed a heuristic 
search with 100 replicates of random taxon 
addition and TBR branch-swapping. Qualitative 
characters excluded from the analysis were poly¬ 
morphic within most populations. Character 
states for the quantitative characters were deter¬ 
mined by searching for gaps within the character 
distribution among specimens that were greater 
than 2 times the average population standard 
deviation (Archie 1985). Most quantitative char¬ 
acters were excluded from the cladistic analysis 
because no character states could be defined. The 
character “locule”, referring to the number of 
locules in the ovary, separates the D. yina complex 
from the outgroup taxa. The morphological data 
matrix is provided in Table 2. 

Cytology 

Chromosome counts were obtained from 
unpublished theses (Weiler 1962; Foster 1972) 
and from numerous specimens deposited at the 
UC and JEPS herbaria as chromosome vouchers 
(Appendix 1). Chromosome number was treated 
as an ordered character. All known chromosome 
counts for D. bacigalupii, D. elegans, D. hicor- 
nuta, D. cuspidata, D. ornatissima and D. 
montana report a single number for each of the 
species (Wood 1961; Foster 1972; Weiler 1962; 
Lammers 1993). All samples of these taxa were 


scored based on chromosome counts reported for 
the species, regardless of whether a count was 
obtained from the population sampled here. The 
only exception is D. hacigainpii sample 585-99, 
which was scored as unknown since the popula¬ 
tion is at the limits of the species range, and no 
chromosome counts were available from the 
vicinity. Chromosome counts for D. concolor 
are // = 8 and n — 9 (Weiler 1962; Lammers 
1993). The samples of D. concolor included here 
fall within the known geographic range of // = 9 
reports for D. concolor (Weiler 1962), and were 
scored as such. Populations of D. yina were 
scored based on the geographic proximity of the 
population to a population with a documented 
chromosome number (indicated in Table 2; 
Weiler 1962; Foster 1972). MacClade version 
3.0 (Maddison and Maddison 1992) was used to 
reconstruct the most parsimonious chromosome 
numbers characterizing each node on trees 
produced from the combined analysis of the 
ITS and 3'P7?K datasets. 

Analyses of Combined Molecular, 
Morphological, and Cytological Data 

A partition-homogeneity test (Farris et al. 
1995) performed in PAUP *4.0 (Swofford 2001) 
confirmed combinability of the ITS plus 3T/7/K 
data (P = 0.247; 1000 replicates, heuristic 
searches with random addition and TBR branch 
swapping), and of the molecular data with the 
morphological and cytological data (P = 0.094). 
Morphological and cytological data were com¬ 
bined as a single partition for the test since 
cytological data consisted of only one character 
(chromosome number). A branch and bound 
search of the combined data was conducted 
under a parsimony criterion. Clade support was 
assessed using 10,000 replicates of the “fast” 
bootstrap option. The morphological data came 
from the same or neighboring populations as the 
sequence data (Table 2). The cytological data 
consisted of chromosome numbers and did not 
include information regarding meiotic configura¬ 
tions of chromosomes in hybrid plants. 

Interfertility 

Information regarding interfertility and cross¬ 
ability among members of the D. yina complex 
comes from Weiler’s unpublished thesis (1962), in 
which he documented the results of numerous 
interspecific crosses within Downingia. His data 
include qualitative assessments of seed set, 
germination, and hybrid condition (e.g., flower¬ 
ing, green, chlorotic, dying in seedling stage), 
some quantitative assessments of pollen stain- 
ability, and analysis of meiotic configurations. 

Information regarding interfertility and cross¬ 
ability within D. yina comes from Foster’s 
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Table 2. Morphological Data Matrix Used for Cladistic Analyses. Sample numbers correspond to 
Appendix 1, with the following prefixes: B = Downingia hacigalupii, E = D. elega?is, Y = D. yimt, M = A niontana, 
C = D. concolor, O = D. ornatissinia, BI = D. bicornuta, CU = D. cuspidata. Characters and states are listed in 
Table 1. The “locule” character refers to the number of locules in the ovary [bilocular (0), unilocular (1)]. For D. 
yina the “chromosome” character refers to the chromosome number based on reports or vouchers from the same or 
a neighboring population, indicated in parentheses. This character was included in the analyses of all data 
combined, but was not included in the analysis of morphological data alone. For some samples, the morphological 
data were combined with the molecular data from a neighboring population, indicated in parentheses. 


Sample 






Character 


6 

8 

10 

13 

16 

17 

Locule 

Chromosome 

B Schn/theis 585-99 

0 

1 

1/0 

0 

0 

1 

1 

? 

B Scimitheis 240-95 

1 

1/0 

1/0 

0 

0 

1 

1 

12 

B Sc/ndtlieis 237-95 

1 

1 

0 

0 

0 

1/2 

1 

12 

B Sc/niltheis 231-95 

1 

1 

1/0 

0 

0 

1 

1 

12 

B Sc/ni/theis 251-95 

1 

1 

0 

0 

1/0 

1/2 

1 

12 

E Schidtheis 243-95 

1 

1 

0 

1 

1 

1 

1 

10 

E Schidtheis 242-95 

0 

1 

0 

1 

1 

1/0 

1 

10 

E Schidtheis 320-96 

0 

1 

0 

1 

1/0 

1 

1 

10 

E Weiter 60138 (Foster 70-15-4) 

0 

1/0 

0 

1 

0 

0 

1 

10 

Y Os'w'cdd & Ahart 3943 

0 

0 

0 

0 

1/0 

1 

1 

? 

Y Schidtheis 247-95 

0 

0 

0 

0 

0 

0 

1 

10 (Foster 70-96-11) 

Y Tracy 3217 

0 

0 

0 

0 

0 

0 

1 

10 (Foster 70-96-11) 

Y. T. Obrien s.n. 

0 

0 

0 

0 

0 

0 

1 

? 

Y Schidtheis 236-95 

0 

0 

0 

0 

0 

0 

1 

12 (Weiler 60207) 

Y D. Barbe 348 

0 

0 

0 

0 

0 

1/0 

1 

12 (Foster, Siskiyou) 

Y Schidtheis 241-95 

0 

0 

1/0 

0 

0 

0 

1 

12 (Foster, Harney) 

Y Schidtheis 584-99 

0 

0 

0 

0 

0 

0 

1 

12 (Foster 70-43-15) 

Y Schidtheis 245-95 

0 

0 

0 

0 

0 

1/0 

1 

12 (Foster 70-43-15) 

Y Weiter 61449 

0 

0 

0 

0 

0 

0 

1 

10 (Weiler 61200) 

Y Schidtheis 581-99 

0 

0 

0 

0 

0 

0 

1 

10 (Weiler 61451) 

Y R. Bacigahipi 7978 

0 

0 

0 

0 

0 

0 

1 

10 (Weiler 61200) 

Y Cook 962 

0 

0 

0 

0 

0 

0 

1 

8 

Y Peck 16291 (Foster 68-210) 

0 

0 

0 

0 

0 

0 

1 

6 (Foster 68-210) 

Y Schidtheis 319-95 

0 

0 

0 

0 

1/0 

0 

1 

12 (Weiler 61383) 

Y Weiter 61333 (Foster 68-51) 

0 

0 

0 

0 

0 

0 

1 

10 

Y R. Bacigahipi 7894 

0 

0 

0 

0 

0 

0 

1 

12 (Foster,Wasco) 

BI Schidtheis 100-95 

0 

0 

1/0 

0 

0 

0 

0 

11 

C Schidtheis 195-95 

0 

0 

0 

0 

0 

1/0 

0 

9 

M Schidtheis 235-95 

0 

0 

0 

0 

0 

0 

1 

11 

CU Schidtheis 179-95 (197-95) 

0 

0 

0 

0 

0 

0 

0 

11 

O Schidtheis 180-95 

0 

0 

0 

0 

0 

1 

0 

12 


unpublished thesis (1972). She documents meiotic 
configurations and pollen stainability for crosses 
between individuals of the same and different 
chromosome races. I assigned each of Foster’s 
parent populations to a molecular cladc, based 
either on sequence data from her voucher 
specimens, or on close proximity of the vouch- 
ered population to a population with sequence 
data. 1 applied an ANOVA to Foster’s raw pollen 
stainability data to examine whether there were 
significant decreases in stainability in hybrids 
between versus within chromosome races, and 
between versus within molecular clades. 

Results 

Cladistic Analyses 

Levels of divergence for the ITS dataset ranged 
from 0.0 to 0.017, excluding outgroups. Analysis 


of ITS data resulted in 104 minimum-length trees 
based on 43 parsimony-informative characters 
(length = 100; Cl = 0.90, 0.83 without uninfor¬ 
mative characters; RI = 0.91). Levels of diver¬ 
gence for the 3'P7iK dataset ranged from 0.0 to 
0.027, excluding outgroups. Analysis of the 
3'tniK dataset resulted in 42 trees based on 19 
parsimony-informative characters (length = 67; 
Cl = 0.94, 0.83 without uninformative charac¬ 
ters; RI = 0.89). Combined molecular analyses 
produced 2 trees based on 43 parsimony-infor¬ 
mative characters (length = 158; Cl = 0.91, 0.77 
without uninformative characters; RI = 0.89). 
Combined molecular, morphological and cyto- 
logical analyses produced 72 trees based on 50 
parsimony-informative characters (length = 176; 
C.I. = 0.87, 0.72 without uninformative charac¬ 
ters; RI = 0.86). 

All analyses (ITS dataset; 3'tmK dataset; 
combined molecular datasets; combined molecu- 
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Fig. 2. The strict consensus of 104 niinimuin-length ITS parsimony trees (length = 100; C.I. = 0.90, 0.83 w/o 
uninformative characters; R.I. = 0.91) produced from a heuristic search with 10,000 replicates of random taxon 
addition and TBR branch swapping. Numbers above the branches indicate bootstrap values generated from 10,000 
replicates of the “fast” bootstrap method. Sample collection numbers correspond to Appendix 1. If only a number 
is indicated, the collection was by Schultheis. Chromosome counts are uniform for all species except Downingia 
yina. Sources for D. yina chromosome counts are indicated in Table 2. 


lar, morphological and cytological datasets) 
except that of morphological data alone resulted 
in three main clades or grades (Figs. 2-5). Clade I 
comprised D. elegans and D. yina pro parte. 
Clade II comprised D. yina pro parte. Clade III 
comprised D. hacigalupii and D. yina pro parte. 
Primary differenees among the trees produced 
from different analyses were the following: (1) 
There was a sister relationship between Clades 1 
and II in trees resulting from analyses of all 
datasets but the ITS dataset, in which Clade II is 
aligned with grade III (Fig. 2). (2) Domiingia 


elegans sample Foster 70-15-4 was resolved as 
part of Clade 1 in all trees except those resulting 
from analysis of the 3'//7/K dataset, in which it 
fell in an unresolved position between Clades I 
and II (Fig. 3). This sample is from Snow 
Mountain, in Lake Co., California, at the 
southern limit of the species range (Fig. 1). (3) 
Downingia hacigalupii sample 585-99 is aligned 
with Clade I in ITS trees (Fig. 2), but is sister to 
other members of Clade III in all other trees. 
Sample 585-99 is from Josephine Co., Oregon, at 
the western periphery of the species range 
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Fig. 3. The strict consensus of 42 minimum-length 3'7/viK intron trees (length = 67; C.I. = 0.94, 0.83 w/o 
uninformative characters; R.I. = 0.89) produced from a heuristic search with 10,000 replicates of random taxon 
addition and TBR branch swapping. Numbers above the branches indicate bootstrap values generated from 10,000 
replicates of the “fast” bootstrap method. Numbers below the branches are decay indices. Sample collection 
numbers correspond to Appendix 1. If only a number is indicated, the collection was by Schultheis. Chromosome 
counts are uniform for all species except Downingia vina. Sources for D. yina chromosome counts are indicated in 
Table 2. 


(Fig. 1). (4) When morphological and cytological 
data are combined with molecular data, the 
resulting trees resolve samples of D. elegans and 
of D. bacigalupii as clades within Clade I and 
Clade III respectively (Fig. 5). Downingia baciga¬ 
lupii sample 585-99, however, is resolved as sister 
to Clade III, and D. elegans sample Foster 70-15- 
4 is unresolved within Clade I. 

The strict consensus of 2556 trees based on six 
parsimony-informative characters (length = 9, 
C.I. = 0.67, RI = 0.88) produced by the cladistic 


analysis of only the morphological data showed 
no resolution (figure not shown). 

Morphometric Analyses 

Downingia yina complex. Univariate analyses 
within the D. yina complex showed that the three 
currently recognized species were significantly 
different from one another for numerous charac¬ 
ters, although ranges overlapped for all charac¬ 
ters (Table 3). Anther angle and the angle of 
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Fig. 4. The strict consensus of two minimum-length trees (length = 158; C.I. = 0.91, 0.77 w/o uninformative 
characters; R.I. = 0.89) from a heuristic search of combined ITS and 3'//7/K intron data, with 10,000 replicates of 
random taxon addition and TBR branch swapping. Numbers above the branches indicate bootstrap values 
generated from 10,000 replicates of the “fast” bootstrap method. Numbers below the branches are decay indices. 
Sample collection numbers correspond to Appendix 1. If only a number is indicated, the collection was by 
Schultheis. Chromosome counts are uniform for all species except Downiugia yina. Sources for D. yina chromosome 
counts are indicated in Table 2. 


divergence between lobes of the lower corolla lip 
in particular distinguished D. elegans and D. 
baeigalupii from D. yina. The former two taxa 
had more sharply bent anthers and less divergent 
lower corolla lobes than D. yina. The filament of 
D. baeigalupii was longer on average than that of 
D. elegans and D. yina. Additionally, D. elegans 
could be distinguished by the lack of yellow 
pignientation on the lower corolla lip. 

PCA analyses of all samples using all data 
showed clear separation among D. elegans, D. 


baeigalupii, and D. yina, particularly when 
principal components 1 and HI were plotted 
(Fig. 6). This separation was also clear when only 
ratio characters were used or when ratio charac¬ 
ters were excluded. Characters of particular 
importance in the PCA analyses were the anther 
angle, the filamcnt/back slit ratio, and the angle 
of divergence between the lower corolla lobes. 
The percent of total variance explained by 
components I, II, and III was 45.2%, 17.5%, 
and 12.0%, respectively. 
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Fig. 5. The strict consensus of 72 minimum-length trees (length = 176; C.I. = 0.87, 0.72 w/o uninformative 
characters; R.l. = 0.86) produced from a branch-and-bound search of combined data. The data matrix included the 
ITS, 3'//7/K, morphological, and cytological data. Numbers above the branches indicate bootstrap values generated 
from 10,000 replicates of the “fast” bootstrap method. Numbers below the branches are decay indices. Sample 
collection numbers correspond to Appendix 1. If only a number is indicated, the collection was by Schultheis. 
Chromosome counts are uniform for all species except Downingia yina. Sources for D. yina chromosome counts are 
indicated in Table 2. 


Cluster analysis (not shown) produced two 
main groups, one with D. yina samples and the 
other with a mixture of D. elegans and D. 
bacigalupii samples. One sample of D. elegans 
(Elilers (& Erlansoii 39) and one sample of D. 
bacigalupii (582-99) together joined at the base of 
the D. yina cluster. 

Variation within Downingia yina. Univariate 
analyses revealed that significant character dif¬ 
ferences were evident between D. yina samples 
from the three main molecular clades, but with 
overlapping ranges (Table 4). No qualitative 
characters could be used to uniquely identify 
the three groups. In general, Clade II samples 
tended to be smaller for most quantitative 
characters measured (Table 4). Samples from 


Clade I tended to have a less sharply bent anther 
and a wider angle of divergence between the lobes 
of the lower corolla lip than did samples from 
Clades II and III. 

PCA analyses of only D. yina samples using all 
data did not clearly distinguish between samples 
from Clades I, II, and III (not shown). Discrim¬ 
inant function analysis of D. yina samples showed 
better separation of the three groups, but with 
areas of overlap (Fig. 7). Characters of particular 
importance in the discriminant function analyses 
were anther length, the angle of divergence 
between the lobes of the lower corolla lip, and 
trichome density on the dorsal anther surface. 

Cluster analysis (not shown) grouped all of the D. 
yina samples together, but did not resolve groups 
corresponding to Clade I, II, and III samples. 
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Table 3. Univariate Statistics for the Downingia yina Complex. Means, standard deviations and ranges 
(in parentheses) are provided for each character within each species. Superscripts indicate groups that are 
significantly different from one another using Tukey multiple comparison tests following ANOVA. Groups with no 
superscript or that share a superscript are not significantly different. 


Character 


D. elegans 



D. yina 

D. hacigalnpii 


(n = 62) 



(n = 317) 

(n = 69) 

Sepal (mm) 

5.31 

± 1.30^ (3.0-8.0) 

4.75 

-F 

1.28” (0.50-10.0) 

5.81 ± 1.80^^ (3.0-10.0) 

Side slit (mm) 

2.33 

± 0.45^ (1.5-3.0) 

4.10 

-F 

0.71” (1.75-6.0) 

2.96 ± 0.54^ (2.0-4.25) 

Back slit (mm) 

4.48 

± 0.75^ (3.0-6.0) 

4.80 


0.89” (2.25-7.5) 

3.87 ± 0.9F (2.25-6.0) 

Upper lobe (mm) 

4.04 

± 1.16^ (2.0-7.0) 

3.71 

-F 

0.93^ (2.0-6.75) 

6.75 ± 1.35” (4.0-11.0) 

Lower lobe (mm) 

6.53 

± 1.82^(3.5-12.0) 

6.18 

± 

1.30^ (3.0-9.5) 

8.40 ± 1.80” (5.0-14.0) 

Filament (mm) 

5.27 

± 1.48^ (2.5-7.5) 

3.17 

-F 

0.92” (1.25-5.75) 

7.48 ± 1.50^' (3.25-9.75) 

Anther (mm) 

2.68 

± 0.46^ (1.5-3.5) 

2.18 

± 

0.36” (1.25-3.0) 

2.89 ± 0.42"' (1.25-3.75) 

Anther angle (degrees) 

84.07 

± 14.40^ (28.0-90.0) 

22.18 


9.11” (0.0-51.0) 

88.80 ± 6.75^ (38.0-90.0) 

Lower angle (degrees) 

9.16 

± 11.92^^ (0.0-50.0) 

53.69 

-F 

12.79” (20.0-90.0) 

12.31 ± 9.88^ (0.0-30.0) 

Horns (mm) 

0.44 

± 0.08^ (0.26-0.75) 

0.41 

-F 

0.1 F (0.13-0.79) 

0.60 ± 0.10” (0.32-0.86) 

Filament/anther 

1.95 

± 0.34^ (1.2-2.55) 

1.45 


0.30” (0.625-2.375) 

2.58 ± 0.35^ (1.6-3.6) 

Side slit/back slit 

0.53 

±0.10^ (0.35-1.0) 

0.86 

-F 

0.09” (0.47-1.3) 

0.79 ±0.13^ (0.5-1.11) 

Filament/back slit 

1.18 

± 0.23^ (0.56-1.75) 

0.66 

-F 

0.18” (0.38-2.0) 

1.97 ± 0.36*^ (1.42-3.56) 

Upper lobe/lower lobe 

0.64 

±0.15^ (0.25-1.0) 

0.61 

-F 

0.14^ (0.31-1.28) 

0.81 ± 0.14” (0.54-1.14) 

Back slit/upper lobe 

1.19 

± 0.35^ (0.6-2.22) 

1.37 


0.41” (0.5-3.11) 

0.59 ±0.18*^ (0.35-1.14) 


Cytology 

Chromosome numbers within the Domiingia 
yina complex appear to correspond to the 
molecular clades identified with ITS and STrnK 
sequences (Figs. 2—4, Appendix 1). All samples in 
molecular Clade I for which chromosome counts 
were available were n — 10 in £>. elegans and n = 
6, 8 or 10 in D. yina. Downingia yina counts of n 


1 -1-r 


"I-i-r 


H 

g 

g 

Cu 

S 

O 

u 

< 

Ph 

O 


ai 

Cu 


2 - 


1 - 


0 - 


1 - 


-2 


a 1=1 




° O o 

pO o 


o 


.O'"© o 

• • 


J_L 


■2-10 1 2 3 4 

PRINCIPAL COMPONENT (3) 


Fig. 6. Plot of principal components one and three 
using the characters listed in Table 1. Filled symbols 
represent Clade I. Symbols with a line through the 
center represent Clade II. Open symbols represent 
Clade III. Triangles = Downingia elegans. Squares = 
D. bacigalupii. Circles = D. yina. 


= 6 and /z = 8 were documented from Marion 
and Lane counties in Oregon, respectively (Weiler 
1962; Foster 1972). All samples in molecular 
Clade II were D. yina with n — 10. All samples in 
molecular Clade III were ;z = 12 in either D. yina 
or D. hacigalnpii. Character state reconstruction 
suggests an ancestral state of zz = 10 in Clades I 
and II, and an ancestral state of zz = 12 in Clade 
III. The ancestral state for the entire D. yina 
complex is equivocal. 

Interfertility 

Foster’s results (1972) show that within D. 
yina, crosses between populations with different 
chromosome numbers showed a significant re¬ 
duction in pollen stainability relative to crosses 
between populations with the same chromosome 
numbers (Table 5; Foster 1972). Similarly, pollen 
stainability was significantly reduced in crosses 
between populations presumed to be from 
different molecular clades relative to those 
presumed to be from the same molecular clades 
(Table 5; Foster 1972). 

Weiler’s results (1962) from interspecific recip¬ 
rocal crosses (Table 6) reflect Foster’s results 
within D. yina in that pollen stainability and 
meiotic irregularities seemed to be affected more 
by differences in chromosome number than by 
species identification {D. elegans, D. hacigalnpii, 
or D. yina). Crosses between D. elegans (n = 10) 
and zz = 10 populations of D. yina, for example, 
produced 10 bivalents and no significant reduc¬ 
tions in pollen stainability (Table 6), in contrast 
to the reduction in pollen stainability for crosses 
within D. yina but between populations of 
different chromosome number (Table 5; Foster 
1972). 
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Table 4. Univariate Statistics Within Downingia yina. Means, standard deviations, and ranges (in 
parentheses) are provided for each character within inferred molecular clades. Superscripts indicate groups that are 
significantly different from one another using Tukey multiple comparison tests following ANOVA. Groups with no 
superscript or that share a superscript are not significantly different. 


Character 


Clade I (n = 100) Clade II (n = 30) Clade III (n = 187) 


Sepal (nmi) 

4.78 

H- 

1.03^® (3.0-7.75) 

4.22 

Side slit (mm) 

4.18 

H- 

0.63^ (2.75-6.0) 

3.40 

Back slit (mm) 

4.86 

-h 

0.78^(3.25-7.0) 

3.74 

Upper lobe (mm) 

3.41 

-h 

0.7 F (2.0-5.75) 

3.44 

Lower lobe (mm) 

6.00 

-h 

0.96"^ (3.25-8.5) 

5.09 

Filament (mm) 

2.99 

-h 

0.62^^ (2.0-4.75) 

2.33 

Anther (mm) 

2.12 

-h 

0.33^^ (1.25-2.75) 

1.75 

Anther angle (degrees) 

18.72 

-h 

8.10” (2.0-40.0) 

22.50 

Lower angle (degrees) 

59.42 

-h 

12.24® (35.0-90.0) 

46.27 

Horns (mm) 

0.38 


0.08® (0.19-0.61) 

0.34 

Filament/anther 

1.42 

-h 

0.23 (1.0-2.0) 

1.37 

Side slit/back slit 

0.86 

± 

0.08^ (0.62-1.07) 

0.91 

Filament/back slit 

0.61 

-h 

0.08® (0.5-0.9) 

0.62 

Upper lobe/lower lobe 

0.58 

-h 

0.12^ (0.33-0.9) 

0.68 

Back slit/upper lobe 

1.50 

-h 

0.44® (0.78-3.11) 

1.17 


0.96® 

(0.5-6.0) 

4.82 

0.49® 

(2.75-4.5) 

4.17 

0.53® 

(3.0-5.0) 

4.94 

0.81® 

(2.0-5.0) 

3.90 

0.99® 

(3.25-6.75) 

6.46 

0.45® 

(1.75-3.5) 

3.40 

0.33® 

(1.25-2.25) 

2.28 

6.ID® (7.0-33.0) 

23.96 

12.40 

^ (20.0-68.0) 

52.06 

0.04® 

(0.26-0.42) 

0.42 

0.30 (0.875-2.2) 

1.48 

0.09® 

(0.77-1.25) 

0.85 

0.08^® (0.47-0.75) 

0.69 

0.13® 

(0.45-0.92) 

0.62 

0.4U 

(0.60-2.22) 

1.34 


1.42^ (2.0-10.0) 
0.73^ (1.75-6.0) 
0.88^ (2.25-7.5) 
0.99^ (2.0-6.75) 
1.39^ (3.0-9.5) 
l.OU (1.25-5.75) 
0.32^ (1.25-3.0) 
9.52^ (0.0-51.0) 
12.08^ (22.0-78.0) 
0.13^ (0.13-0.77) 
0.34 (0.625-2.375) 
0.10^ (0.47-1.31) 
0.22^ (0.38-2.0) 
0.15^ (0.31-1.3) 
0.38^ (0.50-2.5) 


Discussion 

The Downingia yina species complex currently 
comprises three species that are readily distin¬ 
guished from one another on the basis of 
morphological characteristics (Weiler 1962; Ayers 
1993; Fig. 6, Table 3). Downingia elegans and D. 
hacigalupii differ from D. yina in that the anthers 
form a sharp angle relative to the filaments, and 
the lower corolla lobes are relatively parallel 
versus divergent in D. yina. The chromosome 
numbers and the yellow patches on the lower 
corolla lobes readily distinguish D.bacigalupii 



canonical factor 1 

Fig. 7. Plot of canonical factors one and two from 
discriminant function analysis of Downingia yina 
samples using the characters listed in Table 1. Filled 
circles = Clade I. Circles with a line through the center 
= Clade II. Open circles = Clade III. 


from D. elegans. As outlined in the introduction, 
previous workers (Peck 1934, 1937; McVaugh 
1941, 1943) recognized that D. yina may represent 
multiple taxa, which were variously named: D. 
yina Applegate, D. yina Applegate var. major 
McVaugh, D. wiUaniettensis Peck, D. pulcherrima 
Peck. Recent molecular analyses lent merit to 
these interpretations, but sampling within the D. 
yina complex was very limited (Schultheis 2001). 
The additional molecular data presented here 
substantiates these patterns, and demonstrates 
that samples of D. yina fall into three separate 
molecular clades, with D. elegans and D. baciga- 
lupii nested within two of these three clades 
(Figs. 2-^). Taken independently, this paraphy- 
letic or polyphyletic pattern with respect to the 
sequence data from either the nuclear or chloro- 
plast genomes (Figs. 2—4) might only represent 
gene rather than organismal phylogenies (Doyle 
1992; Knox 1998). High resolution molecular 
data are expected to reveal patterns in which 
paraphyletic progenitor species (with respect to 
the molecular data) give rise to monophyletic : 
derivative species (Rieseberg and Brouillet 1994; i 
Graybeal 1995; Olmstead 1995), in this case D. | 

Table 5. Mean Percent Stainable pollen in 
Crosses Between Populations of Downingia ; 
YINA. Significant differences occur for populations i 
with the same versus different chromosome numbers | 
and for populations from the same versus different j 
molecular clades. (Raw data taken from Foster (1972) i 
and reanalyzed). j 


Cross type 

Mean 

SE 

n 

P value 

Chromosome numbers same 

93.4 

8.4 

5 

0.001 

Chromosome numbers differ 

49.9 

5.6 

11 


Same molecular clade 

83.0 

8.3 

7 

0.007 

Different molecular clades 

48.3 

7.3 

9 
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Table 6. Mean Percent Stainability and 
Meiotic Configurations in Interspecific 
Crosses Within the Downingia yina Complex. 
Data taken from Weiler (1962). 



D. elegans 

D. hacigalupii 

D. elegans 

>95% 


n = 10 

D. bacigainpii 

30.8-78.3% 

>95% 

n — \2 

lch3 + 9II + 11 


D. yina 

51.3-79.5% 

>95% 

n = \2 

1211 


D. yina 

>95% 

50-70% 

n = 10 

lOII 

lch3 + 9II + 11 


yina independently giving rise to D. elegans and 
D. hacigalupii. If D. yina populations were 
integrated through gene flow with one another, 
but to the exclusion of D. elegans and D. 
bacigainpii, D. yina would eventually proceed to 
monophyly with respect to the molecular data, 
and the currently recognized species would be 
appropriate, or could be accommodated with 
terms indicating their unresolved or transitional 
status (“metaspecies”, Donoghue 1985; “ferre- 
species”, Graybeal 1995; “plesiospecies”, Olm- 
stead 1995). What is compelling in this example is 
the correspondence of gene geneologies from 
more than one gene with geographic, cytological, 
and interbreeding data; a correspondence that 
makes a case for multiple organismal lineages 
(Avise 1994), and thus multiple species (de 
Queiroz 1998, 1999) within D. yina. 

Geography 

In the D. yina complex, cytological races and 
molecular clades appear to be roughly segregated 
along the Cascade Ranges (Fig. 1). When a 
concordant pattern emerges between phylogenet¬ 
ic and geographic subdivisions of a group, this 
often indicates little to no gene flow among 
subdivisions. This point has been emphasized in 
phylogeographic studies (Avise et al. 1987) and 
has received confirmation from population ge¬ 
netic models (Slatkin 1989). The correspondence 
between molecular clades within the D. yina 
complex and the distribution of these clades to 
either the west or east of the Cascade Ranges is 
striking (Figs. 1-5), and suggests that the moun¬ 
tain range serves as a geographic barrier to gene 
flow. The Cascade Ranges have been recognized 
as a geographic barrier in other contexts, clearly 
affecting differences in climate (Peck 1941; Orr 
and Orr 1996), and floristic composition (Peck 
1941) to the west versus the east. The Klamath- 
Siskiyou region at the California-Oregon border 
is where the striking segregation of D. yina 
molecular clades to the east and west of the 
Cascade Ranges is much less evident (Fig. 1). 


Clade I is found to the west of the Cascade 
Ranges, except that D. elegans extends eastward 
into eastern Washington and Idaho. Clade II is 
localized to a region in the Cascade Range of 
southern Oregon (Fig. 1), in the vicinity of Lake 
of the Woods and Upper Klamath Lake, Oregon, 
and cannot readily be designated as “east” or 
“west”. Molecular clade III is primarily east of 
the Cascades, but extends west into the Klamath- 
Siskiyou region. It is possible that the Klamath- 
Siskiyou region was the source from which the D. 
yina complex dispersed northward to the east and 
west of the Cascades, a scenario similar to 
hypotheses of post-glaciation dispersal presented 
by Whittaker (1961) and Soltis et al. (1997). 

Cytology 

Cytological variation within the D. yina 
complex mirrors the molecular phylogeny and 
the geography for the group, with /; = 12 samples 
primarily east of the Cascade Ranges, and ;/ = 10 
samples primarily to the west. Populations of D. 
yina within Clade I have « = 6 or 8 in the 
northwestern reaches of the species range, an 
observation which prompted Foster (1972) to 
suggest a trend of decreasing chromosome 
numbers as one progressed from the southeast 
to the northwest of D. yina's range. Foster’s 
(1972) proposed explanation for this trend, based 
on meiotic configurations in numerous hybrids 
between the different chromosome races of D. 
yina, was that the races arose through Robertso¬ 
nian translocations producing either a dysploid 
series of reductions from a starting point of n = 
12, or a series of reductions from n = \\ with 
an increase to n = 12. Foster’s work (1972) 
unfortunately did not include D. elegans and D. 
bacigainpii, perhaps because the potential deriva¬ 
tion of these taxa from within D. yina was not 
reflected in the taxonomy. If D. elegans and D. 
bacigainpii arose from within D. yina, as suggest¬ 
ed by the molecular data, the simplest explana¬ 
tion is that they arose from = 10 and /? = 12 
populations of D. yina, respectively. The homol¬ 
ogy of D. elegans and n — 10 D. yina genomes, 
and of D. bacigainpii and n = 12 genomes is 
supported by interfertility data discussed below. 

Interfertility 

If D. yina contains the multiple divergent 
lineages suggested by the molecular data, one 
might expect levels of interfertility to correspond 
with the molecular clades. Indeed, levels of 
interfertility appear to correspond more with 
the molecular clades and the chromosome 
numbers of the populations examined than with 
species identification. For example, individuals of 
D. yina from Clade I show greater interfertility 
with D. elegans than with individuals of D. yina 
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from Clade III (Tables 6 and 7). Similarly, 
individuals of D. yiiui from Clade III show 
greater interfertility with D. hacigahipii than with 
individuals of D. yina from Clades I or II. In sum, 
patterns of interfertility do not appear to 
correspond to the species currently recognized, 
but do appear to correspond to chromosome 
races and molecular data, both of which corre¬ 
spond to geography. 

While levels of fertility may be reduced in 
crosses between chromosome races or molecular 
clades, reproductive barriers are not complete. 
Nor are reproductive barriers complete among 
the three species currently recognized. Popula¬ 
tions exist with hybrids between D. bacigahipii 
and D. yina, and between D. elegans and D. yina 
(Weiler 1962; Schultheis personal observation). 
These populations may either resemble a hybrid 
swarm, with a wide variety of hybrid fonns, or 
may contain readily distinguishable parental 
forms and only a few hybrids (Weiler 1962; 
Schultheis personal observation). Regardless of 
whether reproductive barriers are complete or 
incomplete, the currently recognized species of 
the D. yina complex do not correspond to 
patterns of interfertility within the group. 

Hypothesized Organismal Lineages Within the 
DoH'ningia yina Complex 

In sum, there appear to be three main lineages 
within the D. yina species complex. Members of 
the first lineage (Clade I) are characterized by 
either a “/). yina"''' or “D. elegans'' morphology, 
and are distributed primarily west of the Cas¬ 
cades, with D. elegans extending eastward into 
eastern Washington and Idaho. “D. yina" indi¬ 
viduals are n = 6, 8, or 10. “D. elegans" 
individuals are n = 10. Within this lineage, the 
"D. elegans" members form a clade, excluding 
sample Foster 70-15-4, from the southern periph¬ 
ery of the "D.elegans" range. The scant support 
for the “D. elegans" clade comes from morpho¬ 
logical characters, some of which are polymor¬ 
phic within populations (Table 2). 

The second hypothesized lineage (Clade II), 
localized in the Lake of the Woods region of the 
Cascades in southern Oregon, is characterized by 
a "D. yina" morphology and n = 10. Support for 
this clade comes entirely from molecular charac¬ 
ters. 

Members of the third hypothesized lineage 
(Clade III) are characterized by either a "D. yina" 
or “D. haeigalnpii" morphology, n = 12, and a 
distribution primarily to the east of the Cascades, 
into southwestern Idaho and western Nevada, 
and extending westward into the Klamath/ 
Siskiyou region of southern Oregon and northern 
California. Within this lineage, the “/). hacigaln- 
pii" samples form a clade to the exclusion of 
sample 585-99, from the western periphery of the 


range. The "D. bacigalupii" clade is supported 
only by morphological characters (Table 2). 

Morphological analyses presented here were 
unable to clearly distinguish among D. yina 
samples falling into different molecular clades 
(Fig. 7; Table 4), which largely correspond to 
variation in D. yina chromosome numbers. 
Similarly, Foster (1972) was unable to find 
morphological differences corresponding to the 
chromosome races within D. yina. The chromo¬ 
some races and the molecular clades within D. 
yina are morphologically cryptic. Further exam¬ 
ination of morphology may reveal differences 
missed thus far, but even in the absence of such 
differences, it is desirable to recognize what are 
hypothesized to be organismal lineages. 

Based on the information currently available, I 
choose to recognize five species, with names 
assigned based on nomenclatural priority and the 
phylogenetic placement of the type specimens: D. 
elegans (Lindley) Torrey, D. bacigalupii Weiler, 
D. yina Applegate, D. willaniettensis Peck, and D. 
pulcherrima Peck. Ideally taxon names, including 
species names, should only be assigned to clades 
(Mishler and Donoghue 1982; Misher and 
Theriot 2000). This strict application of a 
phylogenetic species concept only applies full 
species status to D. elegans, D. bacigahipii, and D. 
yina sensu stricto. Downingia willaniettensis and 
D. piilcherrinui comprise the “D. yina" samples 
from Clades I and III respectively. These samples 
are not resolved as clades, but may still be named 
as metaspecies (Donoghue 1985), plesiospecies 
(Olmstead 1995) or ferrespecies (Graybeal 1995). 
An alternative to recognizing five species is to 
recognize a single species, D. elegans (based on 
nomenclatural priority), and five varieties. While 
both alternatives recognize the same taxa, differ¬ 
ing only in the rank applied (species or variety), the 
recognition of five species more clearly emphasizes 
the molecular, cytological and fertility diversity 
within this complex group. In this case, names are 
also available at the species rank whereas new 
names or combinations would be needed if the 
taxa were recognized at the varietal rank. 

Features of the five taxa are summarized in 
Table 7. It is unfortunate that the three species 
previously referred to D. yina (D. yina s.s., D. 
n'illainettensis, and D. piilcherrinui) are morpho¬ 
logically indistinguishable given current informa¬ 
tion. Even those features of most importance in 
discriminate function analysis (included in Ta¬ 
ble 7) show such overlap as to be of minimal use 
for field identification. Weiler (1962) did note 
that D. yina tended to be decumbent in the west 
and erect in the east (which would correspond to 
D. n’illainettensis and D. piilclierrinia respective¬ 
ly), but this can be difficult to detect on 
herbarium sheets. This feature, as well as corolla 
coloration (particularly useful for distinguishing 
D. elegans and D. bacigalupii), is worth noting in 




Table 7. Summary of the Features Distinguishing the Five Species Formerly Classified as Downingia elegans, D. bacigalupii and D. yina. 
* Indicates features identified in discriminant function analysis within D. yina. Mean ± standard deviation; based on measurements from samples in this study 
(Appendix 1). 
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Fig. 8. Map illustrating the distribution of Dowuingia yiua sensn strictn relative to adjacent D. pulcherrinia 
populations in the southern Cascade Range of Oregon. The map does not illustrate D. bacigciliipii, which is also 
found in the pictured region. Triangles = D. yiua s.s. {ii = 10, Clade II). Circles = D. ptdcherruna {n ~ 12, Clade 
III). Downingia yimi s.s. samples are located within the southern tip of the Cascades ecoregion of Oregon (following 
Thorson et al. 2003). 


future collections. Unless reliable features are 
identified, we must rely on geographic location 
for field identification, ideally with confirmation 
from molecular and/or cytological data. At 
present I recommend that specimens collected 
west of the Cascades in Oregon and Washington, 
and west of the North Coast Ranges in California 
are best assigned to D. mllaiuettensis. Specimens 
collected east of the Cascades in Oregon or 
Washington are best assigned to D. pulcherrinia. 
Dowuingia pulcherrinia is also located in the 
Klamath and Siskiyou regions of northern 
California (documented in this study as far west 
as Coffee Creek, just west of Clair Eagle Lake, 
Trinity Co.) and southern Oregon (documented 
in this study as far west as Medford, Jackson 
Co.). Downingia pulcherrinia and D. willanietten- 
sis are generally above and below elevations of 
250 m respectively. Downingia yina seiisii strictn is 
localized to the southern tip of the Cascade 
Range in Oregon. This study documents popula¬ 
tions from the northwestern edge of Upper 
Klamath Lake to Lake of the Woods (Klamath 
Co.). Based on my current understanding of the 
distribution for D. yina, I recommend assigning 
to this taxon any collections found in the 
Cascades ecoregion of southern Oregon (eco¬ 
region as delimited in Thorson et al. 2003), while 
assigning those found in neighboring areas 
outside of this ecoregion to D. pulcherrinia. 
Figure 8 provides a map delimiting the distribu¬ 
tion of D. yina relative to D. pulcherrinia. 


Priorities for refining our current understand¬ 
ing of this species complex include obtaining 
molecular data from additional populations 
(particularly at the limits of species ranges, 
including Washington state) additional sampling 
of cytological variation, and exploration of 
morphological or ecological features to distin¬ 
guish D. yina sensii strictn, D. willaniettensis, and 
D. pulcherrinia. 

Key to Taxa of the Downingia yina 
Species Complex 

la. Anthers abruptly bent, >70° to filaments; 
lower corolla lip lobes ± parallel. 

2a. Corolla 3-colored (blue, white, yellow); 
lower corolla lobes obtuse, mucronate 

. D. bacigalupii 

2b. Corolla 2-colored (blue, white); lower 

corolla lobes acute. D. elegans 

lb. Anthers not or ± bent, <45° to filaments; 
lower corolla lip lobes divergent, not parallel. 

3a. Plants generally east of Cascades, extend¬ 
ing into Klamath Ranges in southern 
Oregon and northern California; generally 
>250 m (but <250 m along Columbia 
River, Washington). 

4a. Localized to southern Oregon Cas¬ 
cades, between northwestern Upper 
Klamath Lake and Lake of the 
Woods, plants at 1200-1510 m. . . D. yina 
4b. East of Cascades, extending into 
Klamath Ranges in southern Oregon 
and northern California, plants gener¬ 
ally at <2000 m. D. pulcherrinia 








2010] 


SCHULTHEIS; CRYPTIC SPECIES IN DOWN IN GlA YIN A 


35 


3b. Plants generally west of Cascades in Oregon 
and Washington, and west of North Coast 
Ranges in California; generally <250 m (but 
650 m on Snow Mountain, Lake Co., 
California). D. wilUnnettensis 
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Sandbergh 287, Douglas Co.: Junction of Crab and Wilson Creeks (UC) 
Eyerdam 1232, Mason Co.: 20 mi SW of Shelton, on shore of small lake (UC) 
Meyer 1551, Thurston Co.: 4 mi W of Olympia along Hwy 101 (UC) 

Sandberg c& Leiberg s.n.. Eastern WA (UC) 
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Schultheis 180-95, Stanislaus Co.: N of Turlock Lake. S side of Barnett Rd, off of Crabtree Rd. (JEPS) AF163360, AF176879 





